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1.   Introduction 

The  "contagious"  binomial  distribution  is  a  discrete  probability 
distribution  which  arises  as  the  steady-state  distribution  of  the 
system  state  variable  in  certain  stochastic  response  models. 
Social  scientists  who  have  used  these  models  have  expressed  a  need 

for  a  limiting  form  of  this  distribution,  but  have  been  unable 

2 
to  derive  it.    The  purpose  of  this  paper  is  to  show  that  the 

desired  limit  for  the  "contagious"  binomial  is  the  beta  distribution. 


In  this  paper  we  first  present  the  general  response  model  which 
has  the  "contagious"  binomial  as  its  steady-state  distribution. 
We  then  develop  the  beta  limit.   Finally,  we  examine  two  major 
types  of  application  of  the  general  model. 

2.   The  General  Model 
In  the  general  formulation  the  basic  model  component  will  be 
referred  to  as  an  element.   The  total  behavior  of  the  element  in 
the  model  will  be  termed  system  behavior.   In  the  discussion  of  the 
applications  of  the  model  presented  below,  the  terms  element  and 
system  will  be  identified  with  the  nomenclature  which  has  been 
used  in  the  application  of  this  model  type  to  particular  stochastic 
response  situations. 
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2. 


Discussion  will  be  facilitated  by  first  considering  the 
notation  which  shall  be  used.   In  the  general  model  an  element 
will  be  uniquely  associated  with  one  of  two  states,  A  or  B.   The 
remainder  of  the  notation  is  as  follows: 
N  -  the  total  number  of  elements  in  the  system. 
1  -  the  number  of  elements  out  of  the  total  of  N  elements  that  are  in 

state  A  at  any  particular  time.  Thus  1  is  a  random  variable  which 

represents  the  state  of  the  system.  Clearly,  0  £  i  i  N. 
O  •  the  inherent  transition  intensity  ,or  propensity  for  an  element 

in  state  B  to  shift  to  state  A. 
p  -  the  inherent  transition  intensity  or  propensity  for  an  element  in 

state  A  to  shift  to  state  B. 
7  -  the  additive  incremental  influence  on  the  transition  intensity  of 

an'  element  toward  the  opposite  state  exerted  by  each  element  in 

the  opposite  state. 

Thus  the  propensity  for  an  element  in  state  B  to  shift  to  state  A  is 

composed  of  two  parts;   an  inherent  propensity,  a,   and  an  additive  attractive 

influence  of  the  i  elements  already  in  state  A,  iy.  Similarly  the  propensity 

for  an  element  to  shift  from  state  A  to  state  B  is  p  +  (N-i)7,  since  there 

i 
are  N-i  elements  currently  associated  with  state  B. 

1 

At  the  level  of  the  elements  the  process  has  two  states.  That  is, 
each  element  is  in  either  state  A  or  state  B.  But  at  the  level  of  the 
system  there  are  N  +  1  states  corresponding  to  the  N  +  1  possible  values  of 
the  system  state  variable,  1.  The  model  gives  rise  to  the  following  system 
of  differential  equations  on  the  system  state  probabilities  at  time  t  which 
are  denoted  by  p.(t): 


dpQ(t)   =  -  N  a  P0(t)   +  [p  +  (N  -l)/]   Px(t)  for  i  =  0  (1) 

~dt 

dPt(t)   m  -    ((N  -   i)    (  a  +  I  7)   +  i   [P  +  (N  -   1)   r]  )  P^*) 
+  (H  -  i  +  1)  [a  +  (i  -  D^lP^Ct) 

+  (i  +  1)    [P  +  (N  -   i   -    1)    7]   Pl+1(t)        for  0<  i<N 

dpn(t)   *  -  N  PpN(t)   +  [a  +  (N  -   1)7]   PN.!(t)  for   i  =  N 

~dt 

The  steady-state  distribution  for  the  state  of  the  system  may  be  found 

by  the  simultaneous  solution  of  the  N  +  1  equations  given  as  (1)  when 

dp^Ct)  =  0  for  i  =  0,  1,  ...  N  and  using  the  fact  that      N 

41—  z  p±(t)-i. 

dt  i=0 


eman  presents  this  steady-state  distribution  as 

■r<">r(H+«-±l>   r<f> 


where         (?  }  =  il  (N-i) ! 

Thus,  at  stochastic  equilibrium  the  distribution  of  i,  the  number  of  elements 

in  state  A,  is  given  by  (2).   The  concern  in  this  paper  is  to  derive  the 

\ 
distribution  of  X  =  i/N  as  N  goes  to  infinity.   That  is,  the  distribution  of 

the  proportion  of  elements  in  state  A  is  sought  as  the  number  of  elements 

goes  to  infinity. 


A. 


3.    The  Beta  Limit 

In  order  to  determine  the  function  which 
will  yield  a  continuous  probability  density  function  for  X,  it  is  first 
necessary  to  consider  what  occurs  as  N  goes  to  infinity  in  such  a  way  that 
X  =  1/N  remains  constant.  The  latter  restriction  is  included  to  ensure 
that  the  cumulative  probability  p[X  <  C],  where  C  is  some  constant  between 
zero  and  one,  remains  constant  as  N  goes  to  infinity. 

For  N  finite ,   the  discrete  mass  function  p(X  =  i/N)  for  X  may  be 
represented  by  a  histogram  such  as  that  presented  in  Figure  1  where 
h(X)  denotes  the  height  of  the  histogram  at  X. 


FIGURE  1 
HISTOGRAM  OP  THE  MASS  FUNCTION  p(X) 
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Note  that  the  width  of  each  interval  is  1/N  and  the  area  of  the  rectangle 
about  X  la  the  probability  that  the  system  has  the  proportion  X  ■  i/N  of 
its  N  elements  associated  with  response  A.  That  is, 

P(X)  »  jjh(x). 

Hence 

h(X)  s  NP(X). 

Suppose  now  that  in  the  interval  (0,1)  more  and  more  elements  are 

packed  in,  that  is,  N  becomes  large.  In  this  case  1/N  will  become  very 

small,  and  in  the  limit  the  h(x)  will  be  eo  close  together  that  they  will 

i  , 
trace  out  a  continuous  curve.  This  continuous  curve  la  the  continuous 

\ 
p.d.f .  of  X  b  i/N  when  N  +*«in  such  a  way  that  X  -  i/N  remains  constant. 

That  is,  where  c  indicates  that  X  ■  1/N  remains  constant, 

Urn  h(x)   -   lim  N   P<X>   p  f(X)  (3) 

N-»  oo  N-»o* 

X=(i/N)=C      X=(i/N)=C 


Thus  (3)  gives  the  function  which  oust  be  taken  to  the  limit  in  order  to 
obtain  the  limiting  probability  density  function  of  X,  f (X) . 

For  any  fixed  N,  there  is  a  direct  equivalence  between  p(X)  in  (3) 
and  the  p  given  by  (2) .  Clearly,  for  some  fixed  N,  the  event  X  =  1/N 
occurs  if,  and  only  if,  the  event  i  occurs.  Hence  for  fixed  N  the  events 
X  =  i/N  and  i  are  equivalent. 

Thus  p(X=i/N)  ■  p..  But  the  steady-state  distribution  of  i 

for  some  fixed  N  has  been  given  in  (2).  Hence  (3)  may  be  written  as 

f  (  X)  »  lira  N  p,  /4) 

N  -+t» 
1/N-C 


,    n/n\  -tc* +  !*<*£-  *>   r(9Lf1) 


lira 

S-T-w  r(a):r(N+aL±i)  ^ 


A  result  on  the  limiting  behavior  of  gamma  functions  when  a  term 
in  the  gamma  argument  increases  without  limit  is  needed  if  the  limit  on  the 

right  hand  side  of  (4)  is  to  be  found.   It  is  known  that  : 


(5) 


lira 

a  ■+  "> 

r(« 
aar 

+  a) 
(a) 

1. 

Furthermore, 

lira 
if 

a  -*•  » 

A(o) 
a(a) 

litn 

a  •+  » 

B(a) 
b(a) 

and  if 

lira 

a  -*■  <* 

a(a)   b( 
'  c(a)    d( 

a) 

k,   ththi 

lira 

A(a)    B(a) 
=  k 

lira   C(a)   lim   D(a)v 
■>  °°   c(a)   a  ■*■  oo  d(a)i 


(6) 


a  •*■  »  C(a)  D(a) 
Using  these  results  in  (4)  we  have 


7. 


£(X)   =  lim  Np  (7) 

MUi=NC 

■ lim    ngLr1)  w*  gg  ^'7  ^h  g*g wP/y  (1  • x)P/y  r(w  [1 : xl) 

nx=i     dp    r(^)  nx    r(NX)  n  (l  -  x)    hn  [I  -  xp  \P'7  n*77    r(N) 

=  T<gL±i)     x<a/7>  -1  d-x)(p/7)  -1 
r(S)  k£) 

But  (7)  Is  Just  the  beta  distribution  with  mean,  a/ (a  4-  0),  and  variance, 

2  7 

ap//(a  +  p)   (a  +  8  a  7) .    Thus  the  infinite  element,  continuous  counter- 
part of  the  "contagious"  binomial  distribution  is  the  beta  distribution  given 
in  (7). 

A.   Applications  of  the  General  Model 

There  have  been  two  principal  classes  of  application  of  the 
general  model.   These  classes  are  considered  below  along  with  an  indication 
of  the  motivation  for  seeking  the  limiting  distribution  in  each  case. 

Class  I.   Models  of  Change  and  Response  Uncertainty 

The  focus  in  this  class  of  applications  has  been  upon  modeling  the 
choice  between  two  alternatives  in  situations  where  a  respondent  chooses 
A  versus  B  with  probability  P(A)  on  any  given  choice  occasion  and  where  P(A) 
may  change  between  successive  choice  occasions.   Thus  the  model  has  beep 
used  in  situations  where  choice  is  stochastically  determined  and  where  the 
probability  of  making  a  particular  choice  is  non-stationary.  This  form  of 

the  general  model  has  been  applied  to  data  from  a  social  psychological 

a 

experiment  and  to  consumer  brand  choice. 


In  these  applications  the  elements  of  the  general  model  are  considered 

to  be  hypothetical  response  elements,  similar  in  nature  to  the  stimulus 

9 
elements  of  stimulus  sampling  theory.   The  system  in  the  general  model  is 

now  the  individual  respondent.   At  any  time  t  an  individual's  probability 

of  making  response  A  versus  response  B  is  given  by  X  =  i/N,  the  proportion 

of  his  N  response  elements  which  are  currently  associated  with  response  A  — 

that  is,  in  state  A.   Thus  in  this  case  (2)  represents  the  steady-state 

distribution  of  an  individual  respondent's  probability  of  making  response  A 

versus  response  B.   In  an  application  to  consumer  brand  choice  response  A 

might  be  a  purchase  of  a  brand  of  particular  interest  to  us  and  response  B 

a  purchase  of  some  other  brand.   Note  that  even  at  stochastic  equilibrium, 

an  individual '8  probability  of  making  response  A  may  change. 

The  desirability  of  letting  the  number  of  response  elements  increase 
without  limit  in  this  case  is  clear.  A  model  which  allows  the  probability 
that  an  individual  will  make  response  A  to  take  on  any  value  between  zero 
and  one  is  more  satisfying  than  one  which  constrains  this  probability  to 
certain  discrete  values.  For  example,  a  model  which  considers  each  individual 
to  have  two  response  elements  implies  that  an  individual's  probability  of 
making  response  k   is  either  0,  1/2,  or  1.  A  model  which  allows  each  j 
respondent  to  have  an  unlimited  number  of  response  elements  allows  this 
probability,  X,  to  be  a  continuous  measure.   Hence,  one  motivation  for 
seeking  the  limit  of  (2)  in  this  case  is  to  develop  the  steady-state 
distribution  for  a  more  theoretically  satisfying  model  of  the  choice 
probabilities. 


In  addition,  this  limiting  result  may  be  of  use  in  estimating  the 
distribution  of  P(A)  across  a  population  of  respondents  all  having  the 
same  a,  g,  and  y.    The  fact  that  the  beta  distribution  is  a  well 
tabled  distribution  makes  this  limiting  result  especially  useful  in  such  cases. 
Class  II  Reward  Models  of  Interpersonal  Influence 

Coleman   has  presented  models  of  the  process  of 

interpersonal  influence  where  the  behavior  of  others  influences  an 

II 

individual  to  exhibit  similar  behavior.   In  these  applications  the  elements 

of  the  general  model  correspond  to  individuals  and  the  system  corresponds  to 
the  group.  The  parameters  a   and  p  correspond  to  the  inherent  propensities 
of  individuals  to  shift  from  activity  0  to  activity  A  and  from  A  to  B, 
respectively.   The  parameter  y   measures  the  influence  which  the  behavior 
of  other  individuals  has  over  the  behavior  of  any  given  individual. 
Coleman  has  utilized  this  form  of  the  model  to  analyze  the  impact  of  inter- 
personal influence  upon  voting  behavior  in  small  groups. 

The  motivation  for  developing  the  limiting  form  of  (2)  in  this  case  is 
one  of  convenience  rather  than  theoretical  desirability.  For  moderately 
large  groups,  say  groups  having  over  twenty-five  members,  (2)  will  have 
N  >  25  and  will  be  computationally  inconvenient.   The  well  tabled  beta  'form 
will  often  prove  to  be  a  useful  approximation  in  such  cases. 

5.   Summary 
\ 
This  paper  has  demonstrated  that  the  beta  distribution  is  the  limit 

distribution  of  the  "contagious"  binomial.   The  structure  of  certain  models 

which  have  the  "contagious"  binomial  as  the  steady-state  distribution  was 

illustrated. 


FOOTNOTES 

1.  Coleman  has  coined  the  term  "contagious  binomial"  to  describe  the 
distribution  considered  in  this  paper. 

2.  Coleman  [3,  p.  413]  expressed  the  need  for  this  limit  but  indicated 
that  he  had  been  unable  to  derive  it. 

3.  For  a  discussion  of  the  notion  of  transition  intensity  in  relation  to 

continuous  time  stochastic  processes  see  [5,  pp. 423-8]. 

4.  See  [3,  pg.  345]. 

5.  This  restriction  is  analogous  to  that  used  in  the  development  of  the 
Poisson  limit  of  the  binomial  distribution.   See  [5]. 

6.  See  [4,  Chapter  9]. 

7.  See  [9,  pg.  216].  ' 

8.  See  [2],  [3,  Chapter  13],  [7,  Chapter  5],  and  [8]. 

9.  For  the  derivation  of  this  model  type  from  an  axiomatic  system 
analogous  to  those  used  in  stimulus  sampling  theory  see  [7,  Chapter  1-3], 

10.   See  [7.  pp.  98-102]  for  a  suggested  use  of  the  limiting  distribution 


in  this  context. 
11.   See  [1]  and  [3,  pp.  343-353]. 
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